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Abstract. Because of their enormous intrinsic brightness blue supergiants are ideal 
stellar objects to be studied spectroscopically as individuals in galaxies far beyond 
the Local Group. Quantitative spectroscopy by means of efficient multi-object spec- 
trographs attached to 8m-class telescopes and modern NLTE model atmosphere tech- 
niques allow us to determine not only intrinsic stellar parameters such as effective 
temperature, surface gravity, chemical composition and absolute magnitude but also 
very accurately interstellar reddening and extinction. This is a significant advantage 
compared to classical distance indicators like Cepheids and RR Lyrae. We describe the 
spectroscopic diagnostics of blue supergiants and introduce two concepts to determine 
absolute magnitudes. The first one (Wind Momentum - Luminosity Relationship) uses 
the correlation between observed stellar wind momentum and luminosity, whereas the 
second one (Flux-weighted Gravity - Luminosity Relationship) relies only on the de- 
termination of effective temperature and surface gravity to yield an accurate estimate 
of absolute magnitude. We discuss the potential of these two methods. 

1 Introduction 

The best established stellar distance indicators, Cepheids and RR Lyrae, suffer 
from two major problems, extinction and metallicity dependence, both of which 
are difficult to determine for these objects with sufficient precision. Thus, in 
order to improve distance determinations in the local universe and to assess the 
influence of systematic errors there is definitely a need for alternative distance 
indicators, which are at least as accurate but are not affected by uncertainties 
arising from extinction or metallicity. It is our conviction that blue supergiants 
are ideal objects for this purpose. The big advantage is the enormous intrinsic 
brightness in visual light, which makes them available for accurate quantitative 
spectroscopic studies even far beyond the Local Group using the new generation 
of 8m-class telescopes and the extremely efficient multi-object spectrographs 
attached to them |2] . Quantitative spectroscopy allows us to determine the stellar 
parameters and thus the intrinsic energy distribution, which can then be used to 
measure reddening and the extinction law. In addition, metallicity can be derived 
from the spectra. We emphasize that a reliable spectroscopic distance indicator 
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Fig. 1. Evolutionary tracks of massive stars in the Hertzsprung- Russell diagram and 
the location of blue supergiants of late B and early A spectral types (shaded box). 
The tracks are from |18| . Solid tracks include the effects of stellar rotation, whereas 
dotted tracks are for non-rotating stars. Solar metallicity has been adopted for the 
calculations. The tracks are labelled by the initial masses (in solar units) at the zero- 
age main sequence (ZAMS). Note that effects of mass-loss due to stellar winds are 
included in the calculations so that actual masses are smaller than ZAMS-masses in 
later stages of the evolution 



will always be superior, since an enormous amount of additional information 
comes for free, as soon as one is able to obtain a reasonable spectrum. 

In this review we concentrate on blue supergiants of spectral types late B to 
early A. These are the the brightest "normal" stars at visual light with absolute 
magnitudes —7.0 > My > —9.5, see (Hj- By "normal" we mean stars evolving 
peacefully without showing signs of eruptions or explosions, which arc difficult 
to handle theoretically and observationally. 

Figure^shows the location of these objects in a Hertzsprung-Russell diagram 
(HRD) with theoretical evolutionary tracks. With initial ZAMS-masses between 
15 and 40 M Q they do not belong to the most massive and the most luminous 
stars in galaxies. O-stars can be significantly more massive and luminous, how- 
ever, because of their high atmospheric temperatures they emit most of their 
radiation in the extreme and far UV. Late B and early A supergiants are cooler 
and because of Wien's law their bolometric corrections are small so that their 
brightness at visual light reaches a maximum value during stellar evolution. 

In the temperature range of late B and early A-supergiants there are also 
always a few objects brighter than My = —9.5 . Generally, those are more exotic 
objects such as Luminous Blue Variables (LBVs) with higher initial masses and 
with spectra characterized by strong emission lines and sometimes in dramatic 
evolutionary phases with outbursts and eruptions. Although their potential as 
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distance indicators is also very promising, we regard the physics of their evolution 
and atmospheres as too complicated at this point and, thus, exclude them from 
our discussion. 

The objects of our interest evolve smoothly from the left to the right in 
the HRD crossing the temperature range of late B and early A-supergiants on 
the order of several 10 3 years ^HJ- During this short evolutionary phase stellar 
winds with mass-loss rates of the order 10 -6 M yr _1 or less p|] do not have 
time enough to reduce the mass of the star significantly so that the mass remains 
constant. In addition, as Fig.^shows, the luminosity stays constant as well. The 
fact that the evolution of these objects can very simply be described by constant 
mass, luminosity and a straightforward mass-luminosity relationship makes them 
a very attractive stellar distance indicator, as we will explain later in this review. 

As evolved objects the blue supergiants are older than their O-star progen- 
itors, with ages between 0.5 to 1.3 x 10 7 years JB]- All galaxies with ongoing 
star formation or bursts of this age will show such a population. Because of their 
age they are spatially less concentrated around their place of birth than O-stars 
and can frequently be found as isolated field stars. This together with their in- 
trinsic brightness makes them less vulnerable as distance indicators against the 
effects of crowding even at larger distances, where less luminous objects such as 
Cepheids and RRLyrae start to have problems. 

With regard to the crowding problem we also note that the short evolutionary 
time of 10 3 years makes it generally very unlikely that an unresolved blend of 
two supergiants with very similar spectral types is observed. On the other hand, 
since we are dealing with spectroscopic distance indicators, any contribution of 
unresolved additional objects of different spectral type is detected immediately, 
as soon as it affects the total magnitude significantly. 

Thus, it is very obvious that blue supergiants seem to be ideal to inves- 
tigate the properties of young populations in galaxies. They can be used to 
study reddening laws and extinction, detailed chemical composition, i.e. not 
only abundance patterns but also gradients of abundance patterns as a function 
of galactocentric distance, the properties of stellar winds as function of chemical 
composition and the evolution of stars in different galactic environment. Most 
importantly, as we will demonstrate below, they are excellent distance indicators. 

It is also very obvious that the use of blue supergiants as tools to understand 
the physics of galaxies and to determine their distances depends very strongly 
on the accuracy of the spectral diagnostic methods which are applied. The at- 
tractiveness of blue supergiants for extragalactic work, namely their outstanding 
intrinsic brightness, has also always posed a tremendous theoretical problem. The 
enormous energy and momentum density contained in their photospheric radia- 
tion field leads to significant departures from Local Thermodynamic Equilibrium 
and to stellar wind outflows driven by radiation. It has long been a problem to 
model non-LTE and radiation driven winds realistically, but significant theoreti- 
cal progress was made during the past decade resulting in powerful spectroscopic 
diagnostic tools which allow to determine the properties of supergiant stars with 
high precision. 
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We describe the status quo of the spectroscopic diagnostics in the follow- 
ing chapters. We will then demonstrate how the spectroscopic information can 
be used to determine distances. We will introduce two completely independent 
theoretical concepts for distance determination methods. The first method, the 
Wind Momentum - Luminosity Relationship (WLR), uses the strengths of the 
radiation driven stellar winds as observed through the diagnostics of H Q as a 
measure of absolute luminosity and, therefore, distance. The second method, 
the Flux- weighted Gravity - Luminosity Relationship (FGLR), determines the 
stellar gravities from all the higher Balmer lines and uses gravity divided by the 
fourth power of effective temperature as a precise measure of absolute magni- 
tude. A short discussion of the potential of these new concepts will conclude the 
paper. 

2 Stellar Atmospheres and Spectral Diagnostics 

The physics of the atmospheres of blue supergiant stars is complex and very 
different from standard stellar atmosphere models. They are dominated by the 
influence of the radiation held, which is characterized by energy densities larger 
than or of the same order as the energy density of atmospheric matter. Another 
important characteristic are the low gravities, which lead to extremely low den- 
sities and an extended atmospheric plasma with very low escape velocity from 
the star. This has two important consequences. First, severe departures from 
Local Thermodynamic Equilibrium (LTE) of the level populations in the entire 
atmosphere are induced, because radiative transitions between ionic energy lev- 
els become much more important than inelastic collisions with free electrons. 
Second, a supersonic hydrodynamic outflow of atmospheric matter is initiated 
by line absorption of photons transferring outwardly directed momentum to the 
atmospheric plasma. This latter mechanism is responsible for the existence of 
the strong stellar winds observed. 

Stellar winds can affect the structure of the outer atmospheric layers sub- 
stantially and change the profiles of strong optical lines such as H Q and EL3 
significantly. The effects of the departures from Local Thermodynamic Equi- 
librium ("NLTE") can also become crucial depending on the atomic properties 
of the ion investigated. A comprehensive and detailed discussion of the basic 
physics behind these effects and the advancement of model atmosphere work for 
blue supergiants is given in |Hj and ^T]. More recent work is described in |29j.|19j 
and [23 • 

For late B and early A-supergiants considerable progress has been made 
during the last four years in the development of a very detailed and accurate 
modelling of the NLTE radiative transfer enabling very precise determinations of 
stellar parameters and chemical abundances, see pQ, [201, [22|, [21], (221 an d |2H- 
Figure El gives an impression of the effort put into the atomic models and corre- 
sponding radiative transfer of individual ions. 
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Fig. 2. Model atoms for Nl (top) and Nil (bottom) as used in the NLTE model 
calculations for late B and early A-supergiants; from |23j 

2.1 Effective Temperature and Gravity 

Effective temperature T c g and gravity logg are the most fundamental atmo- 
spheric parameters. They are usually determined by fitting simultaneously two 
sets of spectral lines, one depending mostly on T e ff and the other on log g. 
Figure indicates how this is done in principle. When fitting the ionization 
equilibria of elements spectral lines of two or more ionization stages have to be 
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Fig. 3. Schematic fit diagram of temperature- and gravity-sensitive spectrai lines in 
the (logp, log T e ff)-plane. Along the dashed curve the computed spectral lines of two 
different ionization stages of one element agree with the observations. Typical ion- 
ization equilibria for late B-supergiants are Si II / ill, Ni/ii, Oi/ii and S 11/ ill, and for 
A-supergiants one can use Mgl/11 and Ni/ii. Along the solid curve the computed pro- 
files of the higher Balmer lines agree with the observations. The intersection determines 
T cff and log g 

brought into simultaneous agreement with observations. At different locations 
in the (logg, log T e ff)-plane this can be achieved only for different elemental 
abundances. Thus, along the fit curve for the ionization equilibrium in Fig. [2] 
the abundance of the corresponding element varies and the intersection with the 
fit curve for the Balmer lines leads to an automatic determination of the abun- 
dance of the element, for which the ionization equilibrium is investigated. (Note 
that the old technique of fitting ratios of equivalent widths of lines in different 
ionization stages and to regard those as being independent of abundance is less 
reliable, since the lines might be on different parts on the curve of growth). 

For A-supergiants the technique has been pioneered by [22] • Most recent ex- 
amples for applications are [23], [22], [21], [221- Examples are given in Figs.0H2l 

The accuracy in the determination of T e fi and logg, which can be achieved 
when using spectra of high S/N and sufficient resolution is astonishing. 
Z\T c ff/T c ff ^0.01 and Z\ log 0.05 are realistic values. 

2.2 Chemical Composition 

The development of very detailed model atoms and using new and very accurate 
atomic data, [Zj, jiiOj . has led to an enormous improvement of the precision to 
which elemental abundances even in extreme blue supergiants can be determined 
H3, [121, E2, [21], 122, E2, 123 • On the average, the uncertainties are now 
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Fig. 4. Fit diagram for the supergiants /?Ori (top) and rjLeo (bottom). The curves 
are parameterized by surface helium abundance y (by number); from |24| 

reduced to 0.1 dex in the abundance relative to hydrogen. Figure displays a 
nice example for the fit of the equivalent widths of CNO lines in blue supergiants. 

The amount of information about chemical elements is impressive. Figure [S] 
gives an overview about the chemical elements the abundances of which can be 
determined from the optical spectra of blue supergiants. Figure [S] shows charac- 
teristic abundance patterns, as they can be derived for supergiants in the Milky 
Way and Fig. displays results for two M31 supergiants. 

2.3 Stellar Wind Properties 

In principle, two types of lines are formed in a stellar wind, P-Cygni profiles with 
a blue absorption trough and a red emission peak and pure emission profiles. 
The difference is caused by the re-emission process after the photon has been 
absorbed within the line transition. If the photon is immediately re-emitted by 
spontaneous emission, then we have the case of line scattering with a source 
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Fig. 5. Temperature and gravity dependence of Mg I (top) and Mg II (bottom) of rj Leo. 
Results from the NLTE computations for the final stellar parameters (full line) are 
compared with synthetic spectra for modified parameters (dotted lines), as indicated, 
against observation (dots); from |24j 

function proportional to the geometrical dilution of the radiation field and a 
P-Cygni profile will result. If the re-emission occurs as a result of a different 
atomic process, for instance after a recombination of an electron into the upper 
level or after a spontaneous decay of a higher level into the upper level or after a 
collision, then the line source function will possibly not dilute and may roughly 
stay constant as a function of radius so that an emission line results. Typical 
examples for P-Cygni profiles are UV resonance transitions connected with the 
ground level, whereas excited lines of an ionization stage into which frequent 
recombination from higher ionization stages occurs will produce emission lines. 
H Q in O-stars and early B-supergiants is a typical example for the latter case. 
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Fig. 6. Temperature and gravity dependence of H 7 of 77 Leo. See Fig. |5] for further 
annotations; from (12 1 j 

However, for late B- and A-supergiants, when Ly a becomes severely optically 
thick and the corresponding transition is in detailed balance, the first excited 
level of hydrogen becomes the effective groundstate and H a starts to behave 
like a resonance line showing also the shape of a P-Cygni profile (for a more 
comprehensive discussion of the line formation process in winds see [H] , |1 1] and 
the most recent review |13)L 

Terminal velocities can be determined very precisely from the blue edges of 
P-Cygni profiles and the red emission wings, normally with an accuracy of 5 to 
10 percent (but see ^Sj for details). In addition, H a profiles normally allow for a 
very accurate (20 percent) determination of mass-loss rates in all cases of O, B, 
and A-supergiants [2Z|, [13 ■ HE], but see and j2H] for details and problems. 

Figure ITU1 gives an impression about the accuracy of the stellar wind spectral 
diagnostics for A-supergiants. 

2.4 Spectral Resolution 

For extragalactic applications beyond the Local Group spectral resolution be- 
comes an issue. The important points are the following. Unlike the case of late 
type stars, crowding and blending of lines is not a severe problem for hot massive 
stars, as long as we restrict our investigation to the visual part of the spectrum. 
In addition, it is important to realize that massive stars have angular momen- 
tum, which leads to usually high rotational velocities. Even for A-supergiants, 
which have already expanded their radius considerably during their evolution 
and, thus, have slowed down their rotation, the observed projected rotational 
velocities are still on the order of 30kms _1 or higher. This means that the in- 
trinsic full half-widths of metal lines are on the order of 1 A. In consequence, 
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Fig. 7. Element abundances derived from individual spectral lines of CNO plotted as 
a function of equivalent width. Top: POii, bottom: r/Leo. Open symbols refer to LTE 
calculations for the line formation, whereas solid symbols show the results of NLTE 
radiative transfer, for neutral (boxes) and ionized species (circles). It is evident that 
LTE fails badly, in particular, for stronger lines. The NLTE results are remarkably 
consistent; from \'2b) 



for the detailed studies of supergiants in the Local Group a resolution of 25,000 
sampling a line with five data points is ideal. This is indeed the resolution, which 
has been applied in most of the work referred to in the previous sections. 

However, as we have found out empirically |24|. degrading the resolution to 
5,000 (FWHM = 1 A) has only a small effect on the accuracy of the diagnostics, 
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Fig. 8. Abundance pattern determined for the two Milky Way supergiants j3 Ori and 
jyLeo, relative to the solar standard |Sj on a logarithmic scale. NLTE [filled symbols) 
and LTE abundances [open symbols) for neutral (boxes), single-ionized (circles) and 
double-ionized (diamonds) species. The symbol size codes the number of spectral lines 
analyzed. Error bars represent la-uncertainties from the line-to-line scatter. The grey 
shaded area marks the deduced stellar metallicity within lcr-errors. The NLTE com- 
putations reveal a striking similarity to the solar abundance distribution, except for 
the light elements which have been affected by mixing with nuclear-processed matter; 
from 



as long as the S/N remains high (i.e. 50 or better). Even for a resolution of 2,500 
(FWHM = 2 A) it is still possible to determine T e ff to an accuracy of 2 percent, 
log g to 0.05 dex and individual element abundances to 0.1 or 0.2 dex. 

References [2], 0, U, [S] and Q2] have used FORS at the VLT with a reso- 
lution of 1,000 (FWHM = 5 A) to study blue supergiants far beyond the Local 
Group. The accuracy in the determination of stellar properties at this rather low 
resolution is still remarkable. The effective temperature is accurate to roughly 4 
percent and the determination of gravity remains unaffected and is still good to 
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Fig. 9. Abundance pattern of the two M 31 A-supergiants 41-3712 and 41-3654. See 
Fig. |S|for further annotations; from 12 II 



0.05 dex (an explanation will be given in Sect. 4). However, at this resolution it 
becomes difficult to determine abundance patterns of individual elements (except 
for emission line stars, see 0]) and, thus, one is restricted to the determination 
of the overall metallicity which is still accurate to 0.2 dex. 

We can conclude that an optimum for extragalactic work is a resolution 
of 2,500 (FWHM = 2 A). Blue supergiants are bright enough even far beyond 
the Local Group to allow the achievement of high S/N with reasonable expo- 
sure times at 8m-class telescopes with efficient multi-object spectrographs at 
this resolution, which then makes it possible to determine stellar parameters 
and chemical composition with sufficient precision. This resolution is also good 
enough to determine stellar wind parameters from the observed H Q profiles, 
since the stellar wind velocities (and therefore the corresponding line widths) 
are larger than 150 km s -1 . 
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Fig. 10. Top: The influence of the mass-loss rate M on the H a profile of the M31 
A-supergiant 41-3654. Two models with M = 1.65 and 2.15 x 10~ 6 Mq yr _1 (dashed- 
dotted, dashed) and otherwise identical parameters are shown superimposed on the 
observed profile. Bottom: The determination of Vaa- Two models with Voa = 200 and 
250 kms -1 (dashed, dashed- dotted) and M adopted to fit the height of the emission peak 
are shown superimposed to the observed profile. All other parameters are identical; 
from |16| 

3 The Wind Momentum— Luminosity Relationship 

The concept of the Wind Momentum-Luminosity Relationship (WLR) has been 
introduced by and [23- It starts from a very simple idea. The winds of 
blue supergiants are initiated and maintained by the absorption of photospheric 
radiation and the photon momentum connected to it. Thus, the mechanical 
momentum flow of a stellar wind Mv^ should be a function of the photon 
momentum rate L j c provided by the stellar photosphere and interior 

M Voo = f{L/c) . (1) 

If this is true and if we are able to find this function /, this would enable 
us to determine directly stellar luminosities from the stellar wind by using the 
inverse relation 

L = r 1 (Mv 00 ) . (2) 
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Fig. 11. Sketch of a blue supergiant irradiating its own stellar wind envelope. L„ is 
the spectral luminosity at frequency v. v is the wind outflow velocity at radius r and 
g is the mass density 



In other words, by measuring the rate of mass-loss and the terminal veloc- 
ity directly from the spectrum we would be able to determine the luminosity 
of a blue supergiant. This is an exciting perspective, because it would give us a 
completely new, purely spectroscopic tool to determine stellar distances. Quanti- 
tative spectroscopy would yield T c g, gravity, abundances, intrinsic colours, red- 
dening, extinction, M and t>oo. With the luminosity from the above relation 
one could then compare with the dereddened apparent magnitude to derive a 
distance. 

In the previous section, we already demonstrated how M and Voo can be 
determined from the spectrum with high precision. Thus, deriving the theoretical 
relationship and confirming and calibrating it observationally will enable us to 
introduce a new distance determination method. An accurate analytical solution 
of the hydrodynamic equations of line driven winds has been provided by 
These solutions were then used by ^0] and |2Z] to exactly derive the relationship. 
Here, we apply a simplified approach, see also jll) . which is not exact but gives 
insight in the underlying physics. 

We assume that the wind is stationary and spherical symmetric and obeys 
the equation of continuity 

M = Anr 2 g(r)v(r) . (3) 

Then we consider the star as a point source of photons irradiating and ac- 
celerating its own stellar wind (see Fig. lTl"|) and calculate the amount of photon 
momentum absorbed by one spectral line 

LL„ 4 (l-e- r «)di' wldth L Vi L v 

c L = ?^ (1 ~ e )dV - (4) 

The first factor on the left side gives the total photon momentum rate provided 
by the star, the second describes the fraction absorbed by one spectral line in an 
outer shell of thickness dr. Tj is the optical thickness of such an outer shell in the 
line transition i. L J , i dt ,wldth is the stellar spectral luminosity at the frequency 
of line i multiplied by the spectral width of the line. This luminosity can in 
principle be absorbed by the line if it is entirely optically thick (i.e. n 3> 1). 
However, depending on the optical thickness only the fraction (1 — e _Ti ) is really 
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absorbed. If we are in the supersonic part of the wind, then the spectral width 
cl^width is not determined by the thermal motion of the ions but rather by the 
increment of the velocity outflow dv via the Doppler formula 

i width 

dv = Vi — , (5) 
c 

which leads to the right hand side of (4). 

After calculation of the photon momentum absorbed by a single spectral line 
we can consider the momentum balance in the stellar wind. The photon momen- 
tum absorbed by all lines will just be a sum over all lines i of the expression 
shown on the right hand side of (4) . Almost all of this absorbed momentum will 
be transformed into gain of mechanical stellar wind flow momentum Mdv of 
the outer shell except the fraction g(r)dM r , which is the momentum required 
to act against the gravitational force (g(r) = GM^/r 2 is the local gravitational 
acceleration and dM r — gAnr 2 dr is the mass within the spherical stellar wind 
shell) 

Mdv = |^^1(1 - e - Ti )dv - G AU(l r ~ r W 2 dr . (6) 

i 

Note that this momentum balance also includes the photon momentum transfer 
by Thomson scattering, which leads to the correction factor (1 — r) in the local 
gravitational acceleration. 

Now, the important next step is to deal with the sum over all lines in the 
above momentum balance. Here, the complication arises from the term in paren- 
theses containing the local optical depth Tj of each line, ti will not only be differ- 
ent for each of the thousands of lines driving the wind, it will also vary through 
the stellar wind as a function of radius. On the other hand, one of the enormous 
simplifications in supersonically expanding envelopes around stars is that the 
optical thickness is well described by (see for instance 8 ) 

, ^thorm ,_n 

Ti = KjK Thom -— — , (7) 

dv/dr 

where Wthcrm is the thermal velocity of the ion and ki is the (dimensionless) line 
strength defined as 



i.e. the opacity of line i 



h = — — , (8) 

KThom 



-mfiu ( 1 - — — ) (9) 



^Dopp m e c \ ni g 



in units of the continuous Thomson scattering opacity or in units of the local 
density 

KThom = n c a c = QS C . (10) 

cr e is the cross section for Thomson scattering of photons on free electrons, n c 
the local number density of free electrons. In a hot plasma with hydrogen as the 
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main constituent mostly ionized and with Yjjc = "He/"H and In e the number of 
electrons provided per helium nucleus we have (mn is the mass of the hydrogen 
atom) 

1 + /He^Hc CT C 

1 + 4r Ho ma 

Thus, if the degree of ionization of helium is roughly constant as a function 
of radius r in the wind, s c is also constant and we have 

Tl = WW 5' (12) 

with the line strength ki proportional to oscillator strength /;„, wavelength A; 
and the occupation number n\ of the lower level divided by the mass density g 

h oc —fi u Xi . (13) 
Q 

Thus, the line strength is roughly independent of the depth in the atmosphere 
and is determined by atomic physics ( fi u , Aj) and atmospheric thermodynamics 
(ni/g). Since g and dv/dr vary strongly through the wind a line can be optically 
thick in deeper layers 

n > 1 =>• (1 - e~ T ' )dv w dv (14) 
and can become optically thin further out 

t, « 1 => (1 - e~ Ti )dv w kigdr . (15) 

For the calculation of the photon momentum transfer this means that line 
contributions can have an entirely different functional form depending on the 
optical thickness of the lines. 

This problem can be solved in a very elegant way by introducing a line 
strength distribution function 

n{k 1 v)dvdk = number of lines with Vi from (y , v + d^) 
and with k from (k, k + dfc) . 

Since the modern hydrodynamic model atmosphere codes (see Sect. 2) con- 
tain atomic data and occupation numbers for millions of lines in NLTE, we can 
investigate the physics of the line strength distribution function. As it turns 
out, PJ, J3| , and the distribution in line strengths - to a very good 
approximation - obeys a power law 

n(k, v)dvdk = g{v)dvk a - 2 dk , 1 < k < oo (16) 

independent (to first order) of the frequency. The exponent a depends weakly 
on T e ff and varies (in the temperature range of OB-stars) between 



a = 0.6... 0.7 



(17) 
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a is mostly determined by the atomic physics and basically reflects the dis- 
tribution function of the oscillator strengths. One can, for instance, show that 
for the hydrogen atom the distribution of the Lyman-series oscillator strengths 
is a power law with exponent a = 2/3, see It is important to realize that 
a is not a free parameter but, instead, is well determined from the thousands of 
lines taken into account in the model atmosphere calculations. Examples for the 
power law dependence of line strengths are given in or . 

With the line strength distribution function we can replace the sum in the 
momentum balance by a double integral, which is then analytically solved (see 
also |H|): 



OO OO 



a-1 



£ ^(1 - e"-) — / /(I - e^f-^n(k, v)dvdk _> N a {^} 

1 o o 

(18) 

This means that the momentum transfer from photons to the stellar wind 
plasma depends non-linearly on the gradient of the velocity field. The degree of 
the non-linearity is determined by the steepness of the line strength distribution 
function a. N a is proportional to the number of lines in the line strength interval 
1 < k < oo. 

We can now re-formulate the momentum balance to obtain 

Mdv = ^N Q dv _ G M,(l-r) 47rrVr (ig) 

yielding a non-linear differential equation for the stellar wind velocity (note that 
we replaced the density through the mass conservation equation) 



^XT = (4-)- 1 r\- - GM*(1 - P) (20) 



2 dv L N ! f 2 dv 

•S— = -. 2-(4tt) q <^ r 2 v— 

dr M a c 2 [ dr 

which looks much more complicated than it really is. The solution is easy, see [3]. 
We obtain M as the uniquely determined eigenvalue of the problem 

MccL^lAf.fl-r)} 1 - 1 /" (21) 

and a terminal velocity proportional to the escape velocity v csc 



Vm OC ZVsc ^X 



{GM t (l-r)/i?,} 1/2 . (22) 



Combining the two yields the stellar wind momentum 



l N L 1 ' a {M,(l-r)} 3 ' 2 - 1 ' a , (23) 



R l/ 



which - as expected - depends strongly on the luminosity but also on the photo- 
spheric radius and the expression in the parentheses, which contains the stellar 
mass and distance from the Eddington limit. It is this expression which can vary 
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Fig. 12. The Wind Momentum - Luminosity Relationship of A-supergiants. Top: ob- 
served stellar wind momentum as a function of absolute magnitude for objects in the 
Milky Way and M31; from Bottom: Theoretical calculations using the stellar wind 
code by |14j . for solar metallicity (solid circles), 0.4 solar (open squares) and 0.1 solar 
(open circles) 



significantly for different blue supergiants and, therefore, causes the large scat- 
ter in the observed correlations of mass-loss rates with luminosity and terminal 
velocity with escape velocity as discussed previously. However, for the product 
of mass-loss rate and terminal velocity, the stellar wind momentum rate, the ex- 
ponent of the term in brackets should be - thanks to the laws of atomic physics 
- very close to zero, since a « 2/3. This means that to first order the wind 
momentum rate should be determined by 

M Voo cx — j-jL 1,/a . (24) 

This is the Wind Momentum - Luminosity Relationship. It predicts a strong 
dependence of wind momentum rate on the stellar luminosity with an exponent 
determined by the statistics of the strengths of the tens of thousands of lines 
driving the wind. It also contains a weak dependence on stellar radius which 
comes from the fact that the stellar wind has to work against the gravitational 
potential when accelerated by photosphcric photons. 
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References [27], and ^2] were the first to prove that the theoretically 
predicted WLR is really observed. O-stars, B and A-supergiants all follow this 
relationship. As to be expected, the relationship depends on spectral types, since 
lines of different ionization stages contribute to the mechanism of line driving at 
different spectra types. F. Bresolin, in this volume [2], gives a detailed overview 
about the most recent observational work on the WLR. Here, we only want 
to show the recent best calibration for A-supergiants using (only) four Milky 
Way and two M31 objects in Fig.^1 Although a calibration based on merely 
six objects is only moderately convincing, we are again encouraged by the small 
scatter over the remarkable range in luminosity. Future calibration work will be 
crucial to establish the method as an accurate distance indicator. 

Also shown in Fig. ^] are recent theoretical stellar wind calculations for 
A-supergiants (R.P. Kudritzki, in prep.), which are based on the new radia- 
tion driven wind algorithm developed by |14j . The calculations provide a clear 
prediction about the metallicity dependence of A-supergiant wind momenta in 
agreement with previous work discussed in and the recent work by |36| 
and F. Bresolin, in this volume |2], will discuss most recent extragalactic 
stellar wind diagnostics on blue supergiants and compare them with the model 
predictions. 

4 The Flux- Weighted Gravity-Luminosity Relationship 

It is an old idea that the strengths of the hydrogen Balmer lines and the abso- 
lute luminosities of massive stars must be related (see 5 in this volume for an 
overview). The concept behind this idea is very simple. Because of the effects of 
Stark broadening the Balmer lines in hot stars are very sensitive to the number 
densities of electrons and protons in those atmospheric layers, where the wings 
of the Balmer lines arc formed. The number densities, on the other hand, depend 
on the photospheric gravity as the result of the hydrostatic equilibrium in stellar 
photospheres. Stellar gravities, however, reflect the evolution of massive stars 
away from the ZAMS (where all gravities are roughly the same) towards lower 
effective temperatures. They become smaller, as further the star evolves, and, if 
we compare objects with exactly the same effective temperature, more massive 
supergiant stars with higher luminosities are expected to have lower gravities 
than their counterparts with lower mass. 

We illustrate the situation by using the stellar evolution models from Fig.^ 
We select three values of T c g — 12000, 9500 and 8350 K, corresponding to spectral 
types B8, AO and A4, respectively (see [2|, this volume), and calculate absolute 
visual magnitudes and gravities for each track at each of the three T e g values. 
Figure H*3l displays the corresponding correlation of My with logg for each tem- 
perature. While the correlations are nicely parallel, the temperature dependence 
as a result of stellar evolution and bolometric correction is obvious. 

We can now calculate atmospheric models for each pair of T e g and log g and 
plot calculated H 7 equivalent widths EW as function of gravity and then, finally, 
produce a diagram of absolute magnitude versus H 7 equivalent width. This is 
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Fig. 13. Three steps to calculate the theoretical correlation between absolute mag- 
nitude My and H 7 equivalent width EW for three different values of T c g = 8350 K 
(crosses), 9500K (squares), 12000K (circles). Top: Mv vs. EW. Middle: EW vs. logg 
as obtained from model atmospheres. Bottom: Mv vs. logg as obtained from stellar 
evolution 



also done in Fig.^J Since the strengths of the Balmcr lines do not only depend 
on gravity but also on temperature, the final correlation 



M v = /(£W(H 7 ) , 



(25) 
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depends very strongly on T c g . The equivalent widths are much stronger at lower 
T e ff and, as an additional complication, the slope of the correlation is strongly 
temperature dependent. While this result is, at least qualitatively, confirmed by 
observation it is clear that simple empirical magnitude - equivalent width 
relations will have to suffer quite some intrinsic scatter unless it is restricted to 
accurate spectral sub-types. We, therefore, suggest a method, which overcomes 
the problem of the strong temperature dependence. 

Before we do this important step, we draw attention to an important detail 
in Fig. The plot equivalent width versus logg reveals that gravities can 
be determined with high precision for each effective temperature. An error of 
^10 percent in EW transforms to an error of 0.05 dex in \ogg. Knowing that 
we will have many higher Balmer lines in the blue spectra of supergiants we 
feel confident that we can determine gravities with such accuracies, even if the 
spectroscopic resolution is only moderate. 

We now turn to the derivation of the Flux- Weighted Gravity - Luminosity 
Relationship (FGLR), which was introduced very recently by [IS]- When dis- 
cussing Fig. Sect. 1 we noted that massive stars evolve through the domain 
of blue supergiants with constant luminosity and constant mass. This has a 
very simple, but very important consequence for the relationship of gravity and 
effective temperature along each evolutionary track. From 

L oc i? 2 T c 4 ff = const.; M = const. (26) 

follows immediately that 

M oc g R 2 oc L (g/T 4 ff ) = const. (27) 

This means that each object of a certain initial mass on the ZAMS has its 
specific value of the "flux-weighted gravity" g/T^ s during the blue supergiant 
stage. This value is determined by the relationship between stellar mass and 
luminosity, which to a good approximation is a power law 

LxM x . (28) 

Inspection of evolutionary calculations with mass-loss, cf. ^7] and US], shows 
that x = 3 is a good value in the range of luminosities considered, although x 
changes towards higher masses. With the mass - luminosity power law we then 
obtain 

L 1 -* oc (g/T^r , (29) 
or with the definition of bolometric magnitude Mboi oc — 2.5 log Z/ 

-M bol = a log( 5 /T c 4 ) + &. (30) 

This is the FGLR of blue supergiants. Note that the proportionality constant 
a is given by the exponent of the mass - luminosity power law through 



a = 2.5x/(l - x) 



(31) 
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Fig. 14. The FGLR of stellar evolution models from |17) and Circles correspond 
to models with rotation, squares represent models without the effects of rotation. Solid 
symbols refer to galactic metallicity and open symbols represent SMC metallicity. The 
solid curve corresponds to a = —3.75 in the FGLR 

and a = —3.75 for x = 3. Mass-loss will depend on metallicity and therefore affect 
the mass - luminosity relation. In addition, stellar rotation through enhanced 
turbulent mixing might be important for this relation. In order to investigate 
these effects we have used the models of 17 and 18 to construct the stellar 
evolution FGLR, which is displayed in Fig. The result is very encouraging. 
All different models with or without rotation and with significantly different 
metallicity form a well defined very narrow FGLR. 

With this nice confirmation of our basic concept we discuss the possible ob- 
servational scatter arising from uncertainties in the determination of T e g and 
log g. As discussed in Sect. 2 and as also demonstrated by |15|. effective temper- 
ature and gravity can be determined within 4 percent and 0.05 dex, respectively. 
Treating these errors as independent we derive an expected one sigma scatter 
zAAlboi = 0.3 mag for the FGLR per individual object, which is again very en- 
couraging and suggests that the method, after careful observational calibration 
(see )15| and jS], this volume), might become a powerful distance indicator. 

Figure ITSI demonstrates how precisely the Balmer lines can be fitted to yield 
a very accurate \ogg and Fig. 1161 shows the first verification of the existence of 
a very tight FGLR for spiral galaxies beyond the Local Group. Further results 
are shown in ^5] and 

5 Conclusions 

We conclude that blue supergiants provide a great potential as excellent extra- 
galactic distance indicators. The quantitative analysis of their spectra - even 
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Fig. 15. Fit of the higher Balmer lines of an A-supergiant in the Sculptor galaxy 
NGC300 using two atmospheric models with T c g — 9500 K and logg = 1.60 (thick 
line) and 1.65 (thin line), respectively. The data were taken with FORS at the VLT. 
For further discussion see \T5\ 

at only moderate resolution - allows the determination of stellar parameters, 
stellar wind properties and chemical composition with remarkable precision. In 
addition, since the spectral analysis yields intrinsic energy distributions over the 
whole spectrum from the UV to the IR, multi-colour photometry can be used to 
determine reddening, extinction laws and extinction. This is a great advantage 
over classical distance indicators, for which only limited photometric information 
is available, when observed outside the Local Group. Spectroscopy also allows 
to deal with the effects of crowding and multiplicity, as blue supergiants, due to 
their enormous brightness, are less affected by such problems than for instance 
Cepheids, which are fainter. 

Two tight relationships exist, the WLR and the FGLR, which can be used 
to derive accurate absolute magnitudes from the spectrum with an accuracy of 
0.3 mag per individual object. Applying the methods on objects brighter than 
My = —8 mag and using multi-object spectrographs at 8 to lOm-class telescopes, 
which allow for quantitative spectroscopy down to my — 22 mag, we estimate 
that with 20 objects per galaxy we will be able to determine distances out to 
distance moduli of m — M ~ 30 mag with an accuracy of 0.1 mag. We note that 
these distances will not be affected by uncertainties in extinction and metallicity, 
because we will be able to derive the corresponding quantities from the spectrum. 
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Fig. 16. The FGLR of B8 to A4 supergiants in NGC 300 and NGC 3621; from P3 
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